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Abstract 
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1 Introduction 



Effective field theory (EFT, for an introduction, see e.g. methods are largely used 
in many branches of physics where a separation of scales exists. In low energy nuclear 
systems, the two well separated scales are, on one side, the low scales of the typical 
momentum of the process considered and the pion mass, and on the other side the higher 
scales associated with chiral symmetry and confinement. This separation of scales was 
explored with great success in the mesonic sector (Chiral Perturbation Theory ^) and 
in the one baryon sector (usually by using Heavy Baryon Chiral Perturbation Theory 
1^), producing a low energy expansion of a variety of observables. It provided for the first 
time a systematic, rigorous and model independent (meaning, independent of assumptions 
about the non-perturbative QCD dynamics) description of strongly interacting particles. 
A large amount of work was devoted recently to the application of EFT methods to 
systems containing two nucleons. In the three nucleon sector, some progress was also made 
in the very low energy range and at low orders where pions do not have to be included 
explicitly f^, 0|. In this article, we investigate the quartet S wave of neutron-deuteron 
scattering using the scheme advanced in ^ in which pions are included perturbatively. 

The original suggestion of how to extend EFT methods to systems containing two or 
more nucleons is due to Weinberg ^ who noticed that below the A production scale, only 
nucleons and pions need to be retained as the infrared relevant degrees of freedom of low 
energy QCD. Because at these scales the momenta of the nucleons are small compared 
to their rest mass, the theory becomes non-relativistic at leading order in the velocity 
expansion, with relativistic corrections systematically included at higher orders. The 
most general chirally invariant Lagrangean consists hence of contact interactions between 
non-relativistic nucleons, and between nucleons and pions. The interactions involving 
pions are severely restricted by chiral invariance. As such, the theory is an extension to 
the many nucleon system of Chiral Perturbation Theory and its counterpart in the one 
nucleon sector. Heavy Baryon Chiral Perturbation Theory. Like in its cousins, all short 
distance physics - quarks and gluons, resonances like the A or a - is integrated out into the 
coefficients of the low energy Lagrangean. In principle, these constants could be derived 
by solving QCD, for example, on the lattice |]T0[. They can also be determined through 
models of the short distance physics, like resonance saturation. The most common way 
to determine those constants, though, is by fitting them to experiment. 

Because its Lagrangean consists of infinitely many terms only restricted by symmetry, 
an effective field theory may at first sight suffer from lack of predictive power. Indeed, as 
part of its formulation, predictive power is ensured only by establishing a power count- 
ing scheme, i.e. a way to determine at which order in a momentum expansion different 
contributions will appear, and keeping only and all the terms up to a given order. The 
dimensionless, small parameter on which the expansion is based is the typical momentum 
Q of the process in units of the scale A at which the theory is expected to break down, 
with estimates ranging from A ^ 300 to 800 MeV [ll|] in the two body system. Assuming 
that all contributions are of natural size, i.e. ordered by powers of Q, the systematic power 
counting ensures that the sum of all terms left out when calculating to a certain order in 
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Q is smaller than the last order retained, allowing for an error estimate of the final result. 

In systems involving two or more nucleons, establishing such a power counting is 
complicated by the fact that unnaturally large scales have to be accommodated, so 
that some coefficients in the Lagrangean may not be of natural size and hence possi- 
bly jeopardise power counting: Given that the typical low energy scale in the problem 
should be the mass of the pion as the lightest particle emerging from QCD, fine tuning 
seems to be required to produce the large scattering lengths in the ^Sq and ^Si chan- 
nels (l/a'^o = -8.3 MeV, l/a^^ = 36 MeV). Since there is a bound state in the ^Si 
channel with a binding energy B = 2.225 MeV and hence a typical binding momentum 
7 = \/MB ~ 46MeV well below the scale A at which the theory should break down, it is 
also clear that at least some processes have to be treated non-perturbatively in order to 
accommodate the deuteron. 

A way to incorporate this fine tuning into the power counting was suggested by Ka- 
plan, Savage and Wise [§]. In it, renormalisation using an unusual subtraction scheme 
(called Power Divergence Subtraction, PDS) moves a somewhat arbitrary amount of the 
short distance contributions from loops to counterterms and makes precise cancellations 
manifest which arise from fine tuning. Power counting then becomes straightforward, 
but physical observables are of course independent of the renormalisation scale or cut-off 
chosen. The leading order contribution to nucleons scattering in an S wave comes from 
four nucleon contact interactions and is summed geometrically to all orders to produce 
the shallow real and virtual bound states. Pion interactions and contact operators con- 
taining derivatives are suppressed by additional powers of and Q, respectively. In the 
pion contributions, the exchange of one instantaneous pion is the leading piece, and pion 
ladders and radiative pions are suppressed further. Dimensional regularisation is chosen 
to explicitly preserve the systematic power counting as well as all symmetries (esp. chiral 
invariance) at each order in every step of the calculation. At leading (LO), next-to- 
leading order (NLO) and often even NNLO in the two nucleon system, this approach also 
allows in general for simple, closed answers whose analytic structure is readily asserted 



1^. This theory has been put to extensive tests at NLO and NNLO in the two body 
system, giving for the first time analytic answers to many deuteron properties such as 
electromagnetic form factors [|1^, scalar and tensor electromagnetic polarisabilities [|14 



Compton scattering [|T^], np — > c/7 both parity violating and conserving [jl^, the 



deuteron anapole moment |T^, charge dependence and charge symmetry breaking 
The expansion parameter is found to be of the order of |, so that NLO calculations are 
accurate to about 10%, and NNLO calculations to about 3%. In all cases, experimental 
agreement is within the estimated theoretical uncertainties, and in some cases, previously 
unknown counterterms could be determined. 

Even if calculations of nuclear properties were possible starting from the underlying 
QCD Lagrangean, effective field theory simplifies the problem considerably by factorising 
it into a short distance part (subsumed into the coefficient of the Lagrangean) and a long 
distance part which contains the infrared-relevant physics and is dealt with by effective 
field theory methods. They provide an answer of finite accuracy because higher order 
corrections are systematically calculable and suppressed in powers of the expansion pa- 
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rameter Q. Hence, the power counting allows for an error estimate of the final result, 
with the natural size of all neglected terms known to be of higher order. The power 
of an effective field theory is also that relativistic effects, chiral dynamics and external 
currents are included systematically, and that extensions to include e.g. parity violating 
effects are straightforward. Gauged interactions and exchange currents are unambiguous. 
Results obtained with EFT are easily dissected for the relative importance of the various 
terms. Because only S'-matrix elements between on-shell states are observables, ambigu- 
ities nesting in "off-shell effects" are absent. On the other hand, because only symmetry 
considerations enter the construction of the Lagrangean, effective field theories are less 
restrictive as no assumption about the underlying QCD dynamics is incorporated. 

This article presents the first systematic effective field theory calculation of three body 
properties with pions. For the theory in which pions are integrated out, calculations 
of three body observables were presented in |0, and S wave quartet neutron-deuteron 
scattering was computed up to NNLO below the deuteron breakup threshold. We extend 
this calculation above threshold and include pions explicitly. 

We choose the quartet S wave channel of nd scattering for this investigation because it 
provides a laboratory in which many complications of the other channels are not encoun- 
tered: The absence of Coulomb interactions ensures that only properties of the strong 
interactions are probed. The Pauli principle forbids three body forces 0] in the first few 
orders. Because the calculation is parameter- free, it allows one to determine the range of 
validity of the KSW scheme without a detailed analysis of the fitting procedure. Although 

the quartet scattering length is large (1/a ^ = 31 MeV), no extra fine tuning except the 
one for the deuteron is required. 

The article is organised as follows: After describing the effective theory and its power 
counting in Sect. |^, we apply it to nd scattering in Sect. |^, followed by a discussion of our 
results in Sect. HI. 



2 The Effective Theory 

2.1 Re-formulating the Lagrangean 

The first terms in the most general Lagrangean satisfying the QCD symmetries including 
pions and nucleons are 

C^^ = iVt(iao + — )iV+ ^tr[(9,St)(9'^S)] + 

+ ^utilMgi^:^ + j:)]+gAN^A-aN+ (2.1) 

- CoiN^P'N)^ (N^P'N) + — \{N^P'Ny{N'^P\d - dfN) + h.c. 

8 L 

- uD2 ti[M^]{N^P'Ny (N^P'N) + ... , 
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where = is the nucleon doublet of two- component spinors and P* is the projector 
onto the iso-scalar-vector channel, 

Paa'' = ^ i^2crr. {r,t , (2.2) 

cr (r) being the Pauli matrices acting in spin (iso-spin) space. The field ^ describes the 
pion, 

e(x) = = e'"/^- , n = 1^ j . (2.3) 

is the chiral covariant derivative = 9^ + V^, and the vector and axial currents are 
given by 

The pion decay constant is normalised to be /^r = 130 MeV, Ai^ = diag(mu,md) is the 
quark mass matrix, and the constant u is chosen such that = u!{^m^^ + m^). From 
now on, we split the coefficient of the leading four point interaction, Cq, into a leading 
(Cq ^■*) and a sub- leading piece (Cg"'*) and absorb D2 into Cg"'* since these two terms can 
be distinguished only in processes involving two or more pions sensitive to explicit chiral 
symmetry breaking by the finite quark masses that will not appear here. 

The Lagrangean has all the symmetries of QCD, including SU{2)i^ x SU{2)^ chiral 
invariance under which the pion and nucleon fields transform as 

^ _ L^U^ = UiR^ , S LTsR^ , N -^UN , (2.5) 

where L and R are constant SU (2) matrices and f/ is a complicated, SU (2) matrix- valued 
function of L, R and the pion fields. The chiral Lagrangean differs from the one used 
in Heavy Baryon Chiral Perturbation Theory by the inclusion of four and other multi- 
nucleon contact forces, the latter of which are not explicitly shown here. 



We now recapitulate the scaling properties of the operators in the Lagrangean ( p.l| ) nec- 
essary to establish the power counting in the next section as suggested by Kaplan, Savage 
and Wise . Because momenta scale like Q and energies like Q'^/M in the non-relativistic 
regime, naive dimensional arguments suggest that the two-nucleon force term with n 
derivatives should scale as ^ 

- ■ (2.6) 

However, the existence of unnaturally shallow bound states and large scattering lengths 
indicates the presence of fine tuning. The inverse deuteron size is 7 = 46 MeV and the 
scattering length in the nucleon-nucleon singlet channel is ~ 24 fm, both of which are 
certainly not typical QCD scales. Most likely, these small scales do not arise from the 
fact that the real world is close to the chiral limit: In the singlet channel, for instance, the 
one pion exchange potential vanishes in the chiral limit and thus cannot be the cause of 
the fine tuning. The fine tuning then must be a result of short distance physics. This, in 
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turn, implies that a naive estimate of sizes of graphs and counterterms using (|2.6| ) may be 
misleading due to precise cancellations and/or enhancements. The exact nature of these 
cancellations and/or enhancements depends on the particular way the short distance 
physics contribution is split between loops and counterterms and, consequently, acquires 
different forms in different regularisation and subtraction schemes. Using a cut-off scheme, 
for example, the sizes of the coefficients C„ satisfy ( p.6|) with A the cut-off scale, but they 
are fine tuned in such a way that sums of diagrams are parametrically larger /smaller than 
each diagram individually, complicating the power counting. It is also hard to keep chiral 
and gauge invariance with those regulators. One convenient way of keeping track of the 
new scale generated by the fine tuning was suggested by Kaplan, Savage and Wise [H. It 
consists in the use of dimensional regularisation and of a subtraction scheme called Power 
Divergence Subtraction (PDS), in which not only the pole present in 4 dimensions is 
subtracted (as in the minimal subtraction scheme), but also poles present in 3 dimensions 
(finite pieces in 4 dimensions) are removed. This unusual subtraction scheme re-shuffles 
a finite piece of the short distance contributions from loops to counterterms and changes 
the scaling of the short distance constants from (|2.6|) to 



(-1) 



MQ 



1 

Co 



MAg2 



Here, Q stands for a low energy scale like the external momenta, the inverse scattering 
length or the pion mass. Other, equivalent ways of accomplishing the same result by 
performing a subtraction at finite off-shell momenta were suggested, too pO| , pT| . One 



surprising result arises from this analysis because chiral symmetry implies a derivative 
coupling of the pion to the nucleon at leading order. The contribution from one pion 
exchange includes a factor of from the pion propagator and a factor of coming 
from the pion-nucleon vertices, so that for momenta of the order of the pion mass, the one 
pion exchange scales as and is smaller than the contact piece Cq"^^ which according 
to ( p.7|) scales as Q~^. Iterated pion exchanges are suppressed even further. On the 
other hand, repeated iterations of the Cq~^^ term are not suppressed, and the two-nucleon 
system is described at leading order by an infinite number of iterations of the leading 
contact interaction, as will be demonstrated below. Pion exchange and higher derivative 
contact terms appear hence only as perturbations at higher orders. In this scheme, the 
only non-perturbative physics responsible for nuclear binding is extremely simple, and 
the more complicated pion contributions are at each order given by a finite number of 
diagrams. In the two body sector, this simplification allows for the derivation of simple 
analytic expressions of deuteron observables. In the three body sector discussed in this 
paper, even the leading order calculation will turn out to be too complex for a fully 
analytical solution. Still, the equations that need to be solved are computationally trivial 
and can furthermore be improved systematically by higher order corrections that involve 
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only (partly analytical, partly numerical) integrations, as opposed to many-dimensional 
integral equations arising in other approaches. 

The estimate of the three-nucleon force terms is also influenced by the fine tuning 
in the two-nucleon sector. In the quartet channel though, where the spin of the three 
participating nucleons are aligned, the Pauli principle forbids a contribution from the 
contact three body force without derivatives. The remaining ones are suppressed and do 
not contribute to the order we are working here. 



For our purposes, it is convenient to use a Lagrangean equivalent to containing an 
additional field carrying the quantum numbers of the deuteron, such that the four 
nucleon contact interactions are removed: 



Nd 



N\ido + 

+ y 



■)N + d'^ 



AM 



w{ido 
+ . . 



where 



4M' 



{2.1 



w 
A(o) 



-1 



(-1) 



MCo 



(2.9) 



(0) 



MCo 



and the pion nucleon interactions are still given by the terms of ( p.l|) . Analogously to 
Co, we split A into a leading (A'^"^-') and a sub-leading piece (A'^^-'). The "wrong" sign in 
the kinetic term for w = — 1 does not spoil unitarity as can be seen from the equivalence 

'-^0 



of (|2.1| ) and (|2.8| ) proven below. Using ( p.9|) and the scaling of the constants Cg ^\ 



and C2 in the KSW scheme, we find 



A(-i) 
A(o) 



y2 ^ 



QA 
M 

Q! 

M 
A 

M2 



(2.10) 



The Lagrangean in ( |2.8|) is not the most general one involving nucleons, pions and 
"deuteron" fields. Since this Lagrangean contains more fields than the minimum number 
necessary, there is a large reparametrisation invariance that can be used to choose the co- 
efficients of some terms to vanish. This freedom in rewriting equivalent Lagrangeans will 
be used later to accomplish non-trivial rearrangements of diagrams in order to simplify 
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computations. To see the equivalence of (|2.8|) and ( pTT|) , we simply perform the GauBian 
integration over the field (P in the path integral. After eliminating the terms containing 
time derivatives by a field re-definition 

^ ^ ^+ 8(AM))2 P'^N\N^P^N) , 
Art ^ Art + ^-^^^ (ArTpw)t™ , (2.11) 



and using ( |2.9|) we are left with ( |2.1|) plus i) terms of higher order, ii) a three-nucleon force 
term with no derivatives that does not contribute to the quartet scattering, and iii) multi- 
nucleon forces not relevant to the three-nucleon system. Because all those contributions 



can be absorbed into terms present already in (|2.1|) the two formulations of the theory are 
therefore indeed equivalent. The choice of ( p.8| ) itself is of course not unique. For example, 
the kinetic term of the deuteron field could have been dropped in favour of keeping the 
C2 term of the nuclear four point interaction. We choose this form for practical reasons 
related to the numerical integration, as will become apparent later. 



2.2 Power Counting and Parameter Determination 

We will briefly state the power counting used here and discuss a couple of examples that 
will be relevant later. More details can be found in [Q. We need to determine the scaling 
of different contributions as the external momenta p, the inverse scattering length a 
and the pion mass go to zero in order to estimate their relative importance. We will 
take p ~ ~ and denote this scale by Q. Our expansion will be in powers of 

Q/K ~ p/A ~ l/(a "^^A) ~ m,r/A, but powers of pa '^^ ~ mj^a ■^i ~ 1 will be kept to all 
orders. 

It is straightforward to estimate the loop contributions: Because the theory is non- 
relativistic, the typical energies of on-shell nucleons are of the order Q^/M. Thus, the 
leading contribution of a diagram can be found by scaling all momenta by a factor of 
Q and all energies by a factor of jM. The remaining integral includes no dimensions 
and is taken to be of the order and of natural size. This scaling implies the rule 
that nucleon propagators contribute one power of M/Q"^ and each loop a power of /M. 
Other regimes, in which loop four momenta scale as Q ~ (on-shell propagation of 
pions in loops) or as ^/Mm^ (pion exchanges with momenta close to the pion production 
threshold), can be shown not to contribute to the order considered here ||2^, |2^. The 
vertices provide powers of Q according to ( p.lO| ), implying that the deuteron kinetic 
energy is sub-leading compared to the A*^^^^ term. Thus, the bare deuteron propagator 
is just the constant —i/A^~^\ Using these rules, we find the diagrams contributing at 
leading order to the deuteron propagator to be an infinite number as shown in Fig. |l], 
each one of the order 1/{MQ). The linear divergence in each of the bubble diagrams 
shown in Fig. |l] does not show in dimensional regularisation as a pole in 4 dimensions, 
but it does appear as a pole in 3 dimensions which we subtract following the PDS scheme 
discussed above. In other regularisation or subtraction schemes, each term in the series 
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Figure 1: The deuteron propagator at leading order from the Lagrangean The 
thick sohd hne denotes the bare propagator the double line its dressed counterpart. 

depicted in Fig. ^ would contribute with a different power of Q, starting with the first, 
tree level diagram that would be of the order 1/(MA). Their sum, however, will still scale 
as 1/{MQ), exemplifying the precise cancellation discussed above, see e.g. ^4 . 



The full leading order propagator i A*-' (p) of the deuteron field consists hence of the 
geometric series shown in Fig. |1|, 

iA'^(p) = - = . (2.12) 

At leading order, all physical quantities will therefore depend on the parameters y and 
A*^~^) only through the combination Cq~^^ = —y'^/^^~^\ so that there is only one free 
parameter in the effective theory. Either y or A'^"^-' can hence be chosen at will. This is 
clear since the same results could have been obtained directly from (p.l|) where only Cq'^^ 



-'0 

appears. We determine this parameter from the binding energy of the deuteron, i.e. from 
the position of the pole of iA(p): 

Eq. (|2.12|) confirms the necessity to re-sum the "bubble chain" with Cq^^^ interactions 
to all orders when considering momenta p ~ 7 and the scaling of Cq'^^ with when 
choosing as subtraction point 11 = ~ Q, in agreement with the claim made previously. 
We will keep /i = throughout this article. By attaching nucleon external lines to 
iA(p), we obtain the nucleon-nucleon scattering amplitude which, at this leading order, 
is simply the effective range expansion truncated at the first term, 

M 7 + 

Looking at ( |2.12| ), we derive the final rule for the power counting: Each deuteron propa- 
gator is accompanied by a power of l/{My'^Q) = M/{AQ). 

At NLO, there are four free parameters in the Lagrangean ( ^.8| ), A'^~^\ A^"), y and 
w, but only three in the equivalent Lagrangean ( p.l|) . Therefore, one of the parameters 
is again arbitrary. We choose to fix the only dimensionless one, w = —1, giving the 
deuteron a kinetic energy term with a supposedly wrong sign. Still, because the two 
Lagrangeans are equivalent, unitarity and causality are not violated. The splitting of 
A = A^"^-* + A''°^ between A*^"^-* and A^") is also arbitrary. We choose to keep A*^~^) the 
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same as at leading order, so that ( p. 13 ) is still valid at NLO and higher order calculations 
will not necessitate a re- fitting of lower order coefficients [^. The remaining piece, 
is then parametrically smaller (as indicated by ( p.lOj )) and is included perturbatively. In 
a NNLO calculation, C2 would be split in a likewise fashion, with 02*^ still given by the 
expression below, C2^^ to be determined anew, and so on. 

The two extra conditions we use to determine the constants A^^^ and y are that the 
position of the deuteron pole 7, which is already reproduced at LO by A*^~^\ is unchanged, 
and that the nucleon-nucleon scattering length in the triplet channel, a is reproduced. 

The condition for the pole position (Fig. |^) implies with w = —1 and fi = (See 
the later, full calculation of the deuteron corrections at NLO in Sect. |3|-) 



A(o) 



11 
M 



(7 



I ln(l 



27 



(2.15) 



327r2/.2 

The first term in brackets comes from the b function part of the pion exchange, the second 
one from the Yukawa part. 




p=0 



Figure 2: The Grst condition on the coefEcients A*^"^-*, A*^'^-' and y, ( |2. J3| / |2!T^ j. The 
dressed deuteron propagator has its pole at the physical binding energy at LO, and higher 
order corrections do not change its position. The cross denotes an energy insertion pro- 
portional to w, the dot an insertion of 

Imposing the experimental value of the triplet nucleon-nucleon scattering length gives 
a second condition (Fig. ^: 

47r AW 



47r 3. 



1 

7 M y'^'-f'^ 



(2.16) 



Pion contributions do not appear in (|2.16|) because they vanish at zero momentum in the 
subtraction scheme used here, the Yukawa piece cancelling against the "5 function at the 
origin" part. Solving ( |2.15| ) and ( p.l6| ) together with ( |2.13| ) leads to 

An 7 

JP G ' 
(m^ - 7)7 J_ 

M G ' 
1 



y 



(2.17) 
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Figure 3: A^A^ scattering at NLO gives the second condition on the coefGcients A^^^ and 



y, ( i2.J6|) ; The scattering length has the experimental value. Notation as in Fig. ^. 



with 



G = 7a - 1 



87r/27 



(7 - m^Y + mi ( ln[l + ^] - 1 



(2.18) 



Figure | shows that the resuhing phase shifts for A^A^ scattering in the ^Si — '^Di system 



in 




200 
k in MeV 





Figure 4: The ^Si and ^Di phase shifts in the NN system and the mixing angle ei 
obtained with the parameters determined as described in the text. The dashed line is the 
LO result, the thin solid (dot-dashed) line the NLO result with perturbative pions (pions 
integrated out), and the thick solid line the Nijmegen phase shift analysis 
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and its mixing parameter is in good agreement with the Nijmegen partial wave anal- 
ysis When pions are integrated out, the phase shift of the ^Di channel and the 
^5*1 — ^Di mixing angle are identical to zero at NLO. At LO, both quantities vanish, too. 



3 Three Body Scattering at NLO 

Let us now look at the diagrams contributing to neutron-deuteron scattering in the quar- 
tet channel. One might guess from the two-nucleon scattering experience that all the 
graphs involving only y and A^"^) (e.g. the ones in the first line in Fig. ^) contribute 
at leading order, and the deuteron kinetic energy, pion exchanges and the A^°^ term ap- 
pear at NLO. Indeed, following the rules discussed above, we can verify that the tree 
level diagram is proportional to y'^M/Q'^ ~ A/(M(5^). The ra-loop diagram with a bare 
deuteron propagator is proportional to ?/2+2n(i/A(-i))"(M/g2)i+2n(Q5/^)n _ a/(^mQ^) 

and hence contributes at the same order. Dressing the deuteron propagator as in Fig. |l] 
does not change the order of the diagram since it only changes the bare propagator 
-i/A(-i) ~ M/(QA) by the dressed propagator iA{p) ~ l/{My^Q) ~ M/(QA). The 
infinite number of graphs shown in Fig. ^ contributing at leading order forms therefore 
a double infinite series. First, all "bubble chain" graphs are re-summed into the leading 




Figure 5: The double inGnite series of LO "pinball" diagrams, some of whicli are shown 
in the first line, is equivalent to the solution of the Faddeev equation ( 3.JS ) shown in the 
second line. Notation as in Fig. 0. 



order deuteron propagator iA(p). Then, the kernel made of iA(p) and the propagator 
of the nucleon exchanged is iterated an arbitrary number of times. We are hence left 
with the task of summing all the diagrams depicted in the first line of Fig. ^. In con- 
tradistinction to the two-nucleon case, they do not form a geometric series and cannot be 
summed analytically. However, one can obtain the solution numerically from the integral 



12 



equation pictorially shown in the lower hne of Fig ^. To derive this equation for the half 
off-shell amplitude at LO, let us choose the kinematics as in Fig. |^ in such a way that 
the incoming (outgoing) deuteron line carries momentum k [pj, energy 3/c^/(4M) —'-f'^/M 
(3fc^/(4M) — + e) and vector index i (j). The incoming (outgoing) nucleon line 

carries momentum —k {—p), energy fc^/(2M) {k^/{2M) — e) and spinor and iso-spinor 
indices a and a {(3 and b). Hence, the incoming deuteron and nucleon are on shell, and e 
denotes how far the outgoing particles are off shell. We denote by it^^^''(/c,p, e) the sum 
of those diagrams with the kinematics above and read off from the lower line of Fig. |^ the 
integral equation for the half off-shell amplitude: 
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(3.19) 
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The integration over go picks the pole at go = (/c^ — (p)/{2M) — e + ie. After that, we set 
e = (/c^ — p^)/(2M), integrate over the angle between k and p to project onto the S wave, 
and set z = (1 + l2)/^/2, j = (1 — i2)/y/2, a = /3 = 1, a = 6 = 1, to pick up the spin 
quartet part. Denoting this projected amplitude by ito(fc,p), we find 
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(3.20) 
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where E = 3k^/{AM) - -/^/M is the total energy. Notice that when p = k, all external 
legs are on-shell. This equation is just the Faddeev equation for the case of contact forces 
derived before by different methods [^. Although the pole in the real axis due to the 
deuteron propagator is regulated by the ie prescription and the logarithmic singularity 
occurring above threshold is integrable, both cause numerical instabilities. We used the 
Hetherington-Schick [^8[ |2^, ^ method to numerically solve (|3.2CI|) . The basic idea is 
to perform a rotation of the variable g into the complex plane by an angle large enough 
in order to avoid the singularities in and near the real axis but small enough so not to 
cross the singularities of the kernel or of the solution. One can show that the singularities 
of the solution are not closer to the real axis than those of the inhomogeneous. Born 
The solution on the real axis can then be obtained from the solution on the 
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deformed contour by using (|3.2CI|) again, now with k and p on the real axis and g on the 
contour. The computational effort then becomes trivial, and a code runs within seconds 
on a personal computer. 



13 



To obtain the neutron-deuteron scattering amplitude, we have to multiply the on-shell 
amplitude tQ{k,p) by the wave function renormalisation constant, 

To{k) - 
1 



Zo to{k, k) sJZq 
d 1_ 



dpQ iA(p) 



Po= 



(3.21) 



In contradistinction to ito{k,p), the scattering amplitude iTo(A;) depends on the parameters 
y and A*^"^^ only through the observable 7. 

The power counting shows that at NLO, we have additional contributions from: deuteron 
kinetic energy insertions, A'-^-' insertions and pions exchange correction to the deuteron 
propagator depicted in the first line of Fig. pionic vertex corrections to y (second and 
third line of Fig. |^); and the pion diagram of the last line of Fig. |^ which corrects the 
three particle intermediate state ("cross diagram"). We call the first three kinds "deuteron 
insertions" and the remaining ones "three-body corrections" : 

iti(fc) = it'^''''''°''{k) + it^^°^y(A;) . (3.22) 



Like the LO deuteron bubbles, the partially diverging off-shell diagrams were calcu- 
lated analytically using dimensional regularisation in order to preserve chiral symmetry 
exactly at each step. The remaining integrations are finite and hence can be treated 
numerically. Each contribution in (|3.22 ) is inserted only once. 



The insertion contributions to the NLO amplitude are given by (see the first line of Fig. | 



it 



insertion 



(k) 



7^ {Mk,q)Y —2 



2M ^ '■^ 



k^ 7^ 



where the corrections to the deuteron are 



lZ^(Po,p) 
with the pion contribution 

(^insertion ^^^^-^ = 
-yinsertion^^^^^^ = 



iAW-i(po-^)+iE.(po,p) 



AM' 



327r2/2 



Y - Mpo - is 



In 




'(Po,p) 



(3.23) 



(3.24) 



(3.25) 
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Figure 6: The NLO contributions to nd scattering in the quartet channel. First hne: 
Corrections to the deuteron propagator; second and third hne: pionic vertex corrections 
to the dNN vertex; fourth hne: pionic corrections to three particle breakup in the inter- 
mediate state (''cross diagrams"). Notation as in Fig. ^. 



Here, as in the following, the pion exchange contributions separate naturally into two 
pieces C^'^^'^^^^°'^{po,p) and 3^™^''''*'°°(po5P) obtained by splitting the pion propagator as 

-»2 2 

1 - ^2 , ^2 • (3-26) 



^2 _|_ j^2 ^2 _|_ ^2 

In configuration space, the first term is a Dirac 6 function resulting in a four nucleon 
contact interaction, and the second is the Yukawa part. As demonstrated at the end of 
this section, this simplifies the calculation further. The 6 function contribution diverges in 
3 dimensions, necessitating a Power Divergence Subtraction of this analytically continued 
pole. The divergence in 4 dimensions of the Yukawa piece has been removed, and the 
arbitrary subtraction constant has been chosen in agreement with p[ such that the total 
pion contribution S^(0,0) at zero momentum and energy is zero. The condition ( p.l5| ) 
that the deuteron pole is not shifted at NLO follows from ( p.24|) , as do the A^A^ phase 
shifts of Fig. |[ 



15 



The go integration in ( p.23| ) picks the nucleon pole, the angular integration is trivial, 
and we are left with a one dimensional integral 



insertion 



dq 
2^ 



iAi 



,3P 
•4M 



21 

M 



1 2 



2M 



(3.27) 



X iX( 



,3^ 
'4M 



11 
M 



2M' 



which has to be performed numerically since we know tQ{k,p) only numerically. 

The computation of the "three-body" diagrams is more involved because a one loop pion 
graph has to be calculated with full off-shell kinematics and then inserted between two 
half off-shell LO amplitudes, resulting after projections in two one dimensional integrals. 
We were able to perform the S wave projection analytically in only part of the graphs 
(the ones in the second and third line of Fig. |]) and we deem it possible to perform also 
the other ones, but the computational advantage will only be minimal since numerical 
integrations have to be performed anyway. Thus, we are left with up to three numerical 
integrals, two over the magnitude of the momenta of the two loops and one over an angle 
in the diagrams of the last line of Fig. |^. 

As depicted in Fig. |^, we choose the four- momentum of the incident deuteron/nucleon 
as {E+e, k)/{^—e, —k), and of the outgoing deuteron/nucleon as {E+e',^))/ (^— e', —p)- 
The parameter e (e') parametrises the off-shellness of the initial (final) state when one 
chooses E = jTf — jTT. The full off-shell behaviour of the one loop pion graphs, projected 



iE+e,k) 




Figure 7: Choice of kinematics for tlie analytical off-sliell calculation of the "three-body' 
diagrams, exemplified at one of the pionic vertex corrections. 

on the quartet S wave state, is then given as 

r(fc,p,e,e',E) = r™^*^-(A;,p,e,e',E) + 



T-'^*-(p, k, + e\ + e,E- ^^) + (3.28) 



2M 

+ r-°-(A;,p,e,e',E) . 



2M 



2M 



The off-shellness of the graphs in the third line of Fig. ^ follows by the kinematic replace- 
ments indicated. 
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The first graph (second hne of Fig. |^) is again spht into a contact and a Yukawa piece 
with the result 



T™^*^^(A;,p,e,e',E) 



— ("P - m^) In 

kp 



r^^'^'^{k,p, e, e', E) - mly^^^'^'^{k,p, e, e', E) 
/C2 + I? + (p - I 



k\2 
2' 



JO + V+{p+\)\ 



(3.29) 



1 

kp 



dz 



1 



In 



+ 2iz/C + ml-K? 
-z"^ + 2izm^ + ml — K? 



- M{E + e)-ie 



where we introduced the convenient abbreviations 

/C 

V 
V 



- M{E + e') - ie 



k^ - p^ 



Me' . 



(3.30) 

(3.31) 
(3.32) 



Again, dimensional regularisation with Power Divergence Subtraction was necessary for 
calculating the contact piece C™''*®^. The analytic answer for the S wave projection rep- 
resented by the integral of the Yukawa piece ( p.30|) is lengthy but straightforward and 
hence not reported here: One uses partial integration followed by a partial fraction of 
the resulting denominator, where the ie prescription takes care that no cuts or poles lie 
in the integration region. The final integrals are standard and lead to logarithms and 
di-logarithms. The S wave projection of the tensor part of the pionic correction to y 
vanishes. 

The contribution corresponding to the full off-shell behaviour of the last line of Fig. ^ 
is finite even in c? = 3: 



r"°^^(A;,p,6,6',i?) 
C"°^^(A;,p,6,e',E) 



1 

kp 



1 

3^^-^^(fc,p,e,e',E) = j 



dy 



4(/C + V) In 
+ Az arctan 

' ^n 



C'°'\k,p, e, 6', E) - m^3;"°-(A;,p, e, e\ E) 
'{k + pf + A{lC + Vf' 



{k-pf + A{lC + Vf 

z 



+ 



kpTZi 
-I- In 



2(/C + V) 
2{k + pf + S 



k+p 
+ 



(3.33) 



2{k-p)^ + S 

2{k + pYyjl -y)+ 47^l7^2(fc + p) + 
2{k - p)^y{l -y)+ 47^l7^2(A; - p) + 
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Here, the following abbreviations are used: 

£ = {2-'iy)k^-{l-?>y)p^ + A{iey + V\l-y))+Am^'R^+Aml{l + y{l-y)) 
T = A{A + y-y')[lC'y + V'{l~y)]+y{l-y)[4ml + 3{k\2-y)+p\l+y) 



7^l = {i2-y)il + y)[2IC''y + 2V''il-y) + mlyil-y) 



+ - y{l-y) [k\2-y)+p\l + y) 
4[IC^y + V\l-y)] + zMl~y) 



+ 



(3.34) 



n2{z) = 

Because it contains polynomials of up to fourth order in y, the Yukawa part y"'°^^ wants 
one numerical integration over the Feynman parameter y - although we deem it possible 
yet time-consuming to find the analytic answer. 

The "three-body" diagrams are classified according to the number of (numerical) con- 
volutions with the LO half off-shell amplitude necessary: 



it?-^°''^(A;) 



.^3-body,o^^^ + 2it?-^°''^''(A;) + it?-^°^^''(A;) (3.35) 



Because the three-body diagrams in the first row of Fig. ^ are not convoluted with the 
LO half off-shell amplitude, 



it 



3- body, 
1 



(k) = ir(fc,fc,0,0. 



4M 



21^ 

M' 



(3.36) 



For the three-body diagrams in the second and third row of Fig. |^ which need to be 
convoluted numerically once and twice with the LO half off-shell amplitude, one picks the 
nucleon pole in each of the integration over the loop energies to obtain 



. ,3— body. 1 
111 



{k) 



dq 
2^ 



ito{k,q) iT{k,q,0. 



k^-q" k' 



2M '4M 



21^ 

M' 



X 



(3.37) 



3-body, 2 



(k) 



X iA 



,3fc2 
'4M 



r 

M 



2M 



(^1 ...2 dp p2 ^^^(^j^^qy^^^(^f.^p-^ X 



27r' 



q 



2tt' 



X iT{p,q,0. 



p^-q^ 3fc2 



2M '4M 



X iAi 



,3^ 



11 
M 



2M 



q) iA 



p^ 

2M ~ 

3k^ 



11^ 
M' 



p 



(3.38) 



AM M 2M 



The wave function renormalisation constant Z at NLO is found from 

1 _ 1 
Z ~ Zo + Zi 
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^0 



1 



P0 = 



r,p = 



as 



Zi - Zq 



(3.39) 
(3.40) 



The NLO amplitude is therefore given by 



T{k) = Zt{k,k) 

- iZo + Z,){toik,k)+hik,k)) 

~ To{k) + Zohik, k) + ZMk. k) = Toik) + T,{k) 



We extract k cot 6{k) from 



37r 

T(A;)^To(A;)+Ti(fc) = — 



1 



M k cot 6{k) - ik 
by expanding k cot 6{k) and keeping only linear terms, i.e. 

1 T^ik)\ 



k cot S 



3m 
ik + — 



M \Toik) Toik)\ 



(3.41) 



(3.42) 



(3.43) 



The phase shift S{k) follows from solving ( p. 431) . Different ways of determining 6 from 
the amplitude yield results which are perturbatively close where the NLO correction 
to the amplitude is parametrically small against the LO amplitude, i.e. below about 
k ^ 200 MeV. 

Finally, we can avoid the computation of the contact pieces of the pion contributions in the 
diagrams of Fig. |^ by exploring the reparametrisation invariance of the Lagrangean ( p. 81) 
in order to get rid of the contact piece of the pion exchange by adding and subtracting a 
term in (p.ll), 



TV AT 



NN 



{N^P'N)\N^P'N) 



A. 



{N'P'Ny{N^P'N) . (3.44) 



The part added is designed to cancel the contact part of the pion exchange and is kept 
in the Lagrangean in ( |2.8| ). The piece subtracted is reproduced in ( p.8| ) by shifting the 
value of the constant A^°) as 



9l 



2/r^i6/2 



{N^P'N)\N^P''N) 



(3.45) 



The NLO wave function renormalisation constant Zi (|3.40|) and the NLO parameters 
(|2.17|/p!T8|) change accordingly. Again, a GauBian integration and field re-definition shows 
that the Lagrangean in ( 3.45| ) is equivalent to the one in ( 3.44|) . 
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At this point, a comment on our choice of the Lagrangean ( p.8|) is in order, too. As 
mentioned above, different Lagrangeans with the fields are equivalent to the Lagrangean 
(|2.1|) . One might, for example, choose to drop the kinetic term of the deuteron field in 
favour of including into ( |2.8| ) the C2 term found in As with the Lagrangean ( |2.8| ), 

the original Lagrangean is obtained when the field is integrated out, showing again 
the equivalence of the three Lagrangeans. In nd scattering, the diagrams to be calculated 
would be different even though their sum produces the same on-shell amplitude. In 
particular, one graph to be calculated is the triangle with C2 on the vertex, dressed 
on one side by the LO solution. Fig. p. According to the analysis of the asymptotic 
behaviour of the solution of the LO integral equation ito{k,p) ~ x/p^'^^ " for p ^ k, 7. 
Consequently, the diagram of Fig. ^ has an UV behaviour of ^q/q2^4 okw = and 
hence is barely convergent, making a numerical integration difficult. In addition, the 
number of diagrams to be calculated is increased considerably. In the approach chosen 




Figure 8: One of the graphs arising in the formulation of the Lagrangean with C2 present, 
having poor convergence properties. 

here, the insertion diagrams have an UV behaviour of Q^^'^^. The three-body diagrams 
which are not convoluted with the LO half of-shell amplitude converge like before, 
and like after absorption of the contact piece of the pion graphs into the definition 
of A^^\ The numerically more involved diagrams of Fig. ^ which have to be convoluted 
once (twice) with the LO half of-shell amplitude converge like Q"^'^^ {Q~^'^^) before, and 
like Q"^-^^ {Q~^-^^) after the contact piece of the pion interaction is removed. 



4 Discussion and Conclusions 

With he = 197.327 MeVfm, a nucleon mass of M = 938.918 MeV, a deuteron binding 
energy (momentum) of 5 = 2.225 MeV (7 = 45.7066 MeV) and the NN triplet S 
wave scattering length a^^^ = 5.42 fm, as well as with the physical values for the pion 
parameters, 171.^ = 138.039 MeV, fj, = 130 MeV and = 1.25, the computation yields 
the results shown in Figs. |10| and |lT] for k cot 6 as function of the momentum in the 
centre-of-mass frame below breakup and for the real and imaginary parts of the phase 



shifts. Because experimental data is scarce in this channel |33] and only available below 



the breakup point /cbrcakup = y^4/3 7 ~ 52.7 MeV, we also compare to the TUNL pd 
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phase shift analysis |^| above breakup as Coulomb and chiral symmetry breaking effects 
can be neglected at high momenta. In addition, two calculations based on the AV 18 
potential @, ^ and the Bonn B potential are presented, the latter one providing 
the only comparison for the imaginary part of the phase shift above breakup. Extending 
the results of Ref. [0] above breakup, we also give an effective field theory calculation to 
NNLO with pions integrated out. Convergence is good. 



breakup 




30 

k in MeV 



Figure 9: k cot S{k) vs. k in the centre-of-mass frame below breakup and extraction of the 
quartet scattering length. The dashed line is the LO result, the solid (dot-dashed) line the 
NLO result with perturbative pions (pions integrated out). The dotted line is an NNLO 
calculation in effective field theory without pions The diamond at zero momentum is 
from the experiment of J3^ . Stars denote the experiment by van Oers and Seagrave as 
reported in |3^/ . Results from realistic potential models are reported as squares from 
crosses from p3 and triangles from 0^ 



In Fig. ^, we focus on the behaviour below threshold where p -C m,T^. In this regime, 
the pions are expected to be integrated out resulting in a much simpler effective ("pion- 
less") theory involving contact interactions between nucleons only 0. Besides our 
NLO result with explicit pions discussed above (full line), we also show the LO, NLO 
and NNLO calculation in a theory with pions integrated out (dashed, dash-dotted and 
dotted line). Simplistically, this is achieved by setting = in the determination of the 
parameters ( p.l7|/^yT^ ) and in the calculations of the NLO corrections. These calculations 
demonstrate convergence and predict a quartet scattering length in excellent agreement 
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breakup 




50 100 150 200 250 300 

k in MeV 

Figure 10: The real part of the phase shift above breakup. The stars are taken from the 
TUNL pd phase shift analysis / |3^/ . Remaining notation as in Fig. ^. 



with experiment. The NLO calculation in the effective theory with and without pions 
give very similar results. The parameters in the pion-less and pion-full LO and NLO 
calculations were determined as described above to produce the physical values of 7 and 
a , while the NNLO result taken from ^ used as input a and the effective range 
To = 1.75 fm, the difference of the two procedures to determine the parameters being of 
higher order in the power counting. No free parameters arise in any of the calculations. 

The difference between the pion-full and pion-less theory should appear for momenta 
of the order of m^r and higher because of non- analytical contributions of the pion cut. 
However, it is very moderate for momenta of up to 300 MeV in the centre-of-mass frame 
{Ecm ~ 70 MeV), see Fig. 110. This and the lack of experimental data makes it difficult to 



assess whether the KSW power counting scheme to include pions as perturbative increases 
the range of validity over the pion-less theory . Unfortunately, numerical calculations with 
realistic potentials are not available for energies considerably higher than the deuteron 
breakup at the present time [|. Still, where data or numerical calculations are available, 
the pion-full theory does better than the pion-less one. 

In deuteron calculations, the numerical value of the systematic expansion parameter Q 
seems to be of the order ^ at zero momentum, and we confirm this by comparing the quar- 



^At the order we are working, an effective theory of a reahstic potential model will be identical to the 
one discussed here. Hence, a model can serve for the present purpose. 
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breakup 




tet scattering length at LO and NLO: We find a^^LO) = (5.1 ± 1.5) fm, and at NLO 

with (without) perturbative pions a 3(NLO,7r) = (6.8 ±0.7) fm (a 5 (NLO, 

(6.7 ±0.7) fm). At NNLO, H report a*^i (NNLO, no tt) = (6.33 ± 0.1) fm, and the 

experimental value is given in [3S] as a ^(exp) = (6.35 ± 0.02) fm. Comparing the NLO 
correction to the LO scattering length provides one with the error estimate quoted: The 
NLO calculation is estimated to be accurate up to about (|)^ ~ 10%. The NLO calcu- 
lations with and without pions lie within each other's error bar. The NNLO calculation 
quoted above is inside the error ascertained to the NLO calculation and carries itself an 
error of about (|)^ ^ 4%. NLO and LO contributions become comparable for momenta 
of more than 200 MeV. In the imaginary part shown in Fig. |TT|, the same pattern emerges 
with a slightly more pronounced difference between the pion-less and pions-fuU theory. 

Because results obtained with EFT are easily dissected for the relative importance 
of the various terms, we show that pionic corrections to nd scattering in the quartet S 
wave channel - although formally NLO - are indeed much weaker: The calculation with 
perturbative pions and with pions integrated out do not differ significantly over a wide 
range of momenta. 



Future work will extend the analysis presented here to higher partial waves and to the 
inclusion of Coulomb effects in order to allow for comparison with pd data. 
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